Sub-Rayleigh Quantum Imaging 
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No imaging apparatus can produce perfect images: spatial resolution is limited by the Rayleigh 
diffraction bound that is a consequence of the imager's finite spatial extent. We show some N- 
photon strategies that permit resolution of details that are smaller than this bound, attaining either 
a 1/y/N enhancement (standard quantum limit) or a 1/N enhancement (Heisenberg limit) over 
standard techniques. In the incoherent imaging regime, the methods presented are loss resistant, 
because they can be implemented with classical-state light sources. Our results may be of importance 
in many applications: microscopy, telescopy, lithography, metrology, etc. 



Quantum effects have been used successfully to provide 
resolution enhancement in imaging procedures. Among 
the many proposals that have been made [1] , arguably the 
most famous is the quantum lithography procedure [2[. 
All of these methods take advantage of the fact that 
the de Broglie wavelength of a multi-photon light state 
is much shorter than the photon's electromagnetic field 
wavelength [3]: the light generation, propagation, and de- 
tection can be performed at optical wavelengths, where 
it is simple to manipulate, whereas the quantum correla- 
tions in the employed states allow one to perform imaging 
at the much shorter de Broglie wavelength. Such propos- 
als are then based on light sources of highly entangled or 
squeezed states, as entanglement or squeezing are neces- 
sary to achieve efficient quantum enhancements [4|]. If, 
however, efficiency considerations are dropped, it is also 
possible to employ classical-state light sources and post- 
selection at the detection stage to filter desirable quan- 
tum states from the classical light [5]. In fact, in many 
practical situations efficiency considerations do not play 
any role, as the quantum enhancement is typically of the 
order of the square root of the number of entangled sys- 
tems [4] , whereas in practical situations the complexity of 
generating the required quantum states has a much worse 
scaling. Many post-selection imaging procedures employ- 
ing only classical light sources have been proposed and 
analyzed @, 0, H !, M, El E H Q EH Efl El H , and 
cover a wide range of interesting situations. Analogous 
methods have been employed successfully also in fields 
not directly related to imaging [19j . 

This Letter discusses how one can achieve a resolution 
enhancement beyond what the apparatus' structural lim- 
its impose for conventional imaging (i.e., the Rayleigh 
diffraction bound xr). In particular we show that em- 
ploying appropriate light sources together with TV-photon 
coincidence photodetection at the output yields a reso- 
lution ~xr/vN. A resolution ~xr/N can also be ob- 
tained by introducing, at the lens plane, a device that 
is opaque when it is illuminated by fewer than N pho- 
tons. The first type of enhancement — a standard quan- 
tum limit for imaging — is an TV-photon quantum process, 
but it is roughly equivalent to the classical procedure of 



averaging the arrival positions of N photons that origi- 
nate from the same point on the object. The second type 
of enhancement — a Heisenberg bound for imaging — is a 
quantum phenomenon that derives from treating the N 
photons as a single field of iV-times higher frequency. In 
the incoherent imaging regime, both methods presented 
here can tolerate arbitrary amounts of loss at the expense 
of reduced efficiency but without sacrificing resolution. 

We start by reviewing some basics of conventional 
imaging. Then we discuss coherent and incoherent sub- 
Rayleigh imaging procedures that attain the standard- 
quantum limit, and finally we introduce our approach to 
realizing the Heisenberg limit for imaging. 

Rayleigh bound:- Consider monochromatic imaging 
using a circular-pupil thin lens of radius R and focal 
length / that is placed at a distance D from an ob- 
ject of surface area *4, and at a distance Di from the 
image plane, where 1/Dq + 1/ Di = 1/ f. In conven- 
tional imaging, the object is illuminated by an appro- 
priate (spatially coherent or incoherent) source and the 
image plane distribution of the light intensity, corre- 
sponding to the probability of detecting a photon at each 
image-plane point f^, is recorded. For photodetectors 
whose spatial-resolution area S and temporal-resolution 
time At are sufficiently small, the preceding probability 
satisfies P^n) ~ (rfScAt)(E^~\fi,t)E^ + \fi,i)), where 
angular brackets denote ensemble average over the illu- 
mination's state, r] is the detector quantum efficiency, 
and E^ — \E^~^ is the positive-frequency compo- 
nent of the electric field. This field component obeys 
E\ + \r u t) = J d 3 k S[ +) \ri-k) e- ikt ^ c a(k), where a(k) is 
the field annihilation operator for the optical mode with 
wave vector fc, and S\ is the solution to the associ- 
ated Helmholtz equation at the image plane. The lat- 
ter can be written in terms of the corresponding object- 
plane field So (r , k) = e lkt ' r ° , where k t is the transverse 
component of fc, using classical imaging equations. For 
monochromatic light in the paraxial regime k t <C /c, it 
follows that (id, El| 



s\ + \n-k) = f^A[f ) h{r u r ) S^\r ;k) 



(1) 



where fi and f are two-dimensional vectors in the im- 
age and object planes, A{f ) is the object aperture func- 
tion [231, and h{fi, f ) is the point-spread function of the 
imaging apparatus given by [13, 120, l2l| 
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with somb(x) = 2J\(x)/x being the well known Airy 
function, and m = Di/D being the image magnifica- 
tion factor. In Eq. (J2J), i? is a phase factor which can be 
generally neglected or compensated. 

Incoherent imaging occurs when the object is illumi- 
nated by independent (monochromatic) beams propagat- 
ing from all directions, whence 
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with I = (Eq Eq ) being the field intensity on the 
object plane. Coherent imaging prevails when collimated 
coherent-state illumination is employed, giving rise to 
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(4) 



When the lens radius R is sufficiently large, Eqs. ([3]) and 
([4]) produce inverted, magnified, perfect images of the ob- 
ject, because R 2 somb(R\x\) — ► 4tt5( 2 \x) for R — > oo [20j . 
For R insufficient to reach this asymptotic behavior, the 
convolution integrals in Eqs. ([3]) and (|4|) produce blurred 
images. The amount of blurring can be gauged through 
the Rayleigh diffraction bound: for a point source at f 
in the object plane, the resulting image-plane intensity is 
proportional to somb (Rk\f + fi/m\/D ), which com- 
prises a pattern of circular fringes in fi that are centered 
on —mf . The radius of the first fringe [211 ] . 



x R ~ 0.61 x 27rmD /(kR) , 



(5) 



about —mf encloses ~84% of the light falling on the im- 
age plane. Intuitively, the image of an extended object 
is then a weighted superposition of radius-x^ circles of 
centered about each —mf . This is the Rayleigh diffrac- 
tion bound; using conventional imaging techniques one 
cannot resolve details smaller than xr. 

Standard quantum limit:- The main idea of sub- 
Rayleigh imaging is to use an appropriate light source 
and to replace intensity measurement with spatially- 
resolving TV-fold coincidence detection strategies. Specif- 
ically we will focus on the probability of detecting TV 
photons at position fi on the image plane [22], i.e., 



N 



(n) - ^^([Et\n,t)] N [E^(n,t)] N ) , (6) 



which can be accomplished by means of doppleron ab- 
sorbers [24|, photon- number resolving detectors, or TV- 
fold coincidence counting. The last two approaches are 
more convenient than the first, as they exploit the full 




FIG. 1: Comparison of conventional, standard quantum limit, 
and Heisenberg limit coherent-imaging for a case in which 
the Rayleigh bound is unable to resolve any object details. 
(a) Object to be imaged, (b) Conventional coherent image 
computed for D /R = 250, m = 1, and k = 6000. The 
length scale for k is given by the image width, and the seg- 
ment in the top right corner is the Rayleigh bound xr. Owing 
to diffraction, no object details are discernible in the image. 
(c) TV-fold coincidence detection (TV = 5) with coherent illu- 
mination, i.e., the reconstruction through Eq. J7}. No res- 
olution enhancement over the previous case occurs despite 
the TVth-order compression of the point-spread function, (d) 
Standard-quantum-limit reconstruction with illumination by 
the superposition of Fock states from Eq. J5J with TV = 5 and 
Akt = 600. Sub- Rayleigh resolution is present, (e) The same 
as (d) except that TV = 10; more resolution enhancement 
occurs, (f) Heisenberg-limited coherent reconstruction from 
Eq. (|12[) with TV = 5; still further enhancement is evident. 



photon statistics so that the TV value need not be prede- 
termined. Note that multi-photon detection alone does 
not guarantee sub- Rayleigh performance. In fact, for the 
coherent imaging of Eq. (|4j) , TV-photon detection gives 
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Here, the factor of TV in the exponent gives an TV- fold 
compression of the fringes in the point-spread function. 
This compression, however, is not intrinsically quantum. 
It amounts to taking the TVth power of the light intensity, 
which is simply a classical post-processing of the signal 
in Eq. ([4]). Thus no resolution enhancement is obtained 
through simple TV-photon detection, see Fig. [T{c). 

As our first example of a source that can be used to 



beat the Rayleigh bound, consider an input state that 
is the superposition of TV-photon Fock states that have 
been focused to a small area s = (nAk 2 ) -1 <C A centered 
at positions r on the object plane, viz., 
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where b(f* ) is the annihilator of the associated localized 
spatial mode [25[ and A4 ~ T^pr is a normalization con- 
stant. Inserting this state into Eq. ([6]), we find 
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Q(n,r ) =f£f A(r) h(n,r) F Akt (\r - r\) , 



where F^kt( x ) 
V 



irAk 2 A somb(Afe* x/2), and £ 



TrAfeM ^ * s a dimensionless quantity that is typically 



very small because of the monochromatic (AcuAt <C 1) 
and focusing assumptions (irAk^A ^> 1). Equation ([9]) 
can be simplified by assuming D /R ^> k/ Ak tl which 
implies that each number state in the superposition is 
focused to a spot much smaller than the object-plane 
Rayleigh limit of the lens. In this case h can be extracted 
from the integral yielding Q(ri,r ) ~ h(fi,f )A(f ) with 
A{f ) = J ^A(r)F Akt (\r - f|). Now Eq. ® becomes 
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which, contrary to Eq. Q, cannot be obtained through 
post-processing of Pi, and which generalizes coherent 
imaging (|4]) to iV-photon detection. The point-spread 
function that governs spatial resolution is now h N — 
which is narrower than h — so that when A(f ) c± A(f ) 
there is an enhancement in resolution over the Rayleigh 
bound. More generally, even if A and A differ signif- 
icantly, one can still beat the Rayleigh bound if N is 
sufficiently large and Di/R ^> k/Ak t: see Figs. QJd) and 

An analogous generalization for incoherent imaging is 
obtained by replacing the state Eq. (|8]) with an incoherent 
mixture of focused Fock states, i.e., p = J ^-^\N)^ o (N\. 
In this case Eq. ([10|) becomes 
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which generalizes Eq. ([3|) to A"-photon detection. The 
corresponding resolution enhancement is shown in Fig.[2l 
The states employed in Eqs. ([8j) and (fTTj) are highly 
sensitive to loss. Nevertheless, A 7 "- fold incoherent imag- 
ing can be realized with loss-resistant light sources. Sup- 
pose we use an incoherent mixture of coherent states 
that randomly illuminate all points on the object: a = 
f ^a 2 l a )r H, where \a)? o =exp[atf(r ) — a*b(r o )]\0). 



Equation ([TT]) still applies with an extra multiplicative 
factor of \a\ 2N /N\ to account for the Poissonian pho- 
todetection statistics of coherent states. The state a can 
be prepared by shining a highly-focused laser beam on 
the object, one point at the time. This state is highly 
robust to loss, because loss parameter fi just takes \a) 
into \y/Jia). Hence an arbitrary amount of loss can be 
tolerated — without sacrificing resolution — simply by in- 
creasing \a\. 
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FIG. 2: Comparison of conventional and standard quantum 
limit incoherent imaging of the object shown in Fig.[TJa) with 
the same D Q /R = 250, m = 1, and k = 6000. (a) and (b) Con- 
ventional images from Eq. (0 with N = 1 and 5, respectively. 
The images are featureless blurs, because the Rayleigh bound 
xr is unable to resolve any object details, (c) and (d) Re- 
construction via Eq. (|11|) using the state a for N = 5 and 
10, respectively, with Ak t = 600. An obvious increase in 
resolution is seen. 

The improved resolution afforded by the procedures 
detailed above can be roughly estimated by gauging by 
narrowing of the point-spread function h that results 
from taking its Nth power. This can be done, for in- 
stance, by evaluating the radius xr(N) that contains 84% 
of the area under somb in the plane. Numerical anal- 
ysis shows that xr(N)/xr ~ l/\/~N, which suggests a 
standard quantum limit [l| for imaging. This should be 
taken only as a rough estimate, as xr(N) is also the ra- 
dial dimension of a point-like object imaged using the 
post-processing strategy of Eq. J7j). For more extended 
objects, the actual resolution enhancement will also de- 
pend on Akt. The 1/y/N scaling exposes the classical 
nature of this enhancement: the same effect can be at- 
tained by averaging the arrival positions of N photons 
at the image plane. This is surely advantageous over N- 
photon detection in many situations, but it is impractical 
for lithography or film photography, and it seems impos- 
sible to classically reproduce the coherent imaging case 
of Eq. (pTU]) . In addition, from general principles [l|,|4| one 
would expect the ultimate bound (a Heisenberg limit for 



imaging) to have 1/A 7 ", scaling, i.e., a resolution ~xr/N, 
which is not achievable with classical strategies. 

Heisenberg limit:- The Heisenberg 1/N scaling can 
be obtained by treating the N photons as a single entity 
of TV-times higher frequency. This situation can be sim- 
ulated, at least in principle, by inserting immediately in 
front of the lens a screen divided into small sections each 
of area sf such that if less than N photons reach one 
section, they are absorbed, otherwise they are coherently 
transmitted. Such a screen does not currently exist, but 
in principle one could be built, e.g., by using doppleron 
materials [24]]. Then, if the object is illuminated by the 
focused coherent states described above, only N photons 
that originate at r , successfully transit the screen within 
one of its area-s^ segments, and get detected at rj can 
can contribute to the image at that point. In this case, 
the operators [E\ \ri,t)} N of the iV-photon absorption 
probability ([TO]) are approximately [26[ 
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Here: Hn is obtained from Eq. ([2]) by replacing k with 
Nk, i.e., hjsf is the point-spread function for photons hav- 
ing TV-times higher frequency than the illumination; and 
7 accounts for the spatial resolution of the doppleron 
screen, i.e., it is of order (^^■) iV - Equation ([T2|) describes 
the absorption of N frequency-cj photons that originated 
near f and then propagated through the imaging ap- 
paratus as if they were a single frequency- Nuo photon. 
It gives rise to coherent and incoherent images that are 
formally equivalent to those of Eqs. (|4j) and (|3]) for a 
light beam of wave number Nk, thus realizing the Heisen- 
berg limit of an TV-fold resolution improvement over the 
Rayleigh bound, see Fig. DJf). This method shows that 
Heisenberg-limited resolution can be obtained using clas- 
sical light, albeit with an even worse efficiency than the 
TV-photon detection methods given above. 
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where the frequency integration is performed over 
a small bandwidth Au around the mean frequency 
uj — kc and the wave-vector integration is performed 
over \kt\ ^ Akt. The exact expression for M is 
fcPfofcPfo F^ro-^/A", which, for highly lo- 
calized modes irAk^A < 1, gives M ~ ^rzr- 
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E? ) {ft\t)---E'? , {?f>,t), with Ef>{n,t) being 
the electric field on the screen. The screen removes 
all contributions save those for which the rf fall in 
the same area-s^ region. Equation ([T2]) follows by 
writing the resulting operator as a convolution of 
£o +) (ri 1} , t) ■ • • #o +) (fi N) , t) through the object-to- lens 
transfer function, and by noticing that for focused 
sources it can be approximated by \b(fo )] • 



